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PREFACE 


The  Agri-Environmental  Indicator  (AEI)  Project  of  Agriculture  and  Agri-Food  Canada 
(AAFC)  was  initiated  in  1993  in  response  to  recommendations  made  by  several  agencies, 
organizations  and  special  studies.  The  overall  objective  of  the  project  is  to  develop  and 
provide  information  to  help  integrate  environmental  considerations  into  decision-making 
processes  of  the  agri-food  sector. 

The  project  aims  to  develop  a  core  set  of  regionally-sensitive  national  indicators  that  build 
on  and  enhance  the  information  base  currently  available  on  environmental  conditions  and 
trends  related  to  primary  agriculture  in  Canada.  The  soil  organic  carbon  component  of  the 
Soil  Degradation  Risk  indicator  is  an  important  part  of  the  agri-environmental  indicator  set. 
Indicators  are  also  being  developed  for  other  aspects  of  soil  degradation  risk  and  in 
relation  to  issues  of  water  quality,  agroecosystem  biodiversity,  farm  resource 
management,  agricultural  greenhouse  gases  and  agricultural  production  efficiency. 

Research  results  in  the  form  of  discussion  papers,  scientific  articles  and  progress  reports 
are  released  as  they  become  available.  A  comprehensive  report  is  planned  for  fiscal-year 
1998-1999  which  will  include  data  from  the  1996  Statistics  Canada  Census  of  Agriculture. 

Comments  and  questions  on  this  paper  should  be  addressed  to: 

Dr.  Greg  Wall 

Ontario  Land  Resource  Unit 

Greenhouse  and  Processing  Crops  Research  Centres 

Research  Branch,  Agriculture  and  Agri-Food  Canada 

Guelph,  Ontario 

N1H3N6 

Telephone:  (519)  826-2086 
Facsimile:  (519)826-2090 
E-mail:  wallg@em.agr.ca 
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Abstract 

In  order  to  develop  a  soil  organic  carbon  (SOC)  indicator  for  Eastern  Canada  the 
Century  model  was  run  for  all  Soil  Landscape  of  Canada  (SLC)  polygons  in  the 
Mixedwood  Plains  Ecozone  of  Ontario.  Within  each  polygon,  SOC  dynamics  were 
estimated  for  three  or  four  crop  rotations  and  two  tillage  practices.    It  was  determined  that 
the  Mixedwood  Plains  Ecozone  was  losing  SOC  at  the  rate  of -39.6,  -34.9,  -32.1,  and  - 
24.3  kg  C  ha1  y1  for  the  years  1981,  1986,  1990,  and  1991,  respectively.  Such  rates  of 
loss  are  not  measurable  on  a  year  to  year  basis  and  are  usually  considered  to  be  at  near 
equilibrium.  The  introduction  of  no-till  agriculture  is  the  principal  factor  that  reduced  the 
rate  of  loss  of  SOC. 

There  was  greater  loss  of  SOC  near  the  Southern  tip  of  Ontario  where  agriculture  is 
most  intensive  (less  hay  and  pasture).  The  greatest  loss  of  SOC  was  found  in  corn  and 
soybean  rotations  and  there  was  a  gain  in  SOC  under  hay  and  pasture  rotations.  It  was 
found  that  soil  texture  influences  SOC  dynamics  with  more  SOC  being  lost  in  coarser 
textured  soils. 

The  Century  model  did  a  reasonable  job  in  estimating  SOC  dynamics  in  Ontario. 
The  model  was  developed  and  parameterized  in  the  western  US  grasslands  so  it  was 
expected  to  be  more  sensitive  to  western  Canada  SOC  dynamics.  Considering  that  the 
model  has  been  tested  successfully  worldwide  under  various  climatic,  soil,  and 
management  practices  there  is  likely  no  better  tool  for  developing  a  soil  carbon  indicator  at 
the  present  time. 


1.0  Introduction 

Soil  organic  matter  is  important  because  it  affects  soil  productivity,  the  physical 
well-being  of  soil,  and  the  atmospheric  content  of  C02  and  other  greenhouse  gases. 
Some  recent  studies  have  indicated  that  the  rate  of  change  of  soil  organic  carbon  (SOC) 
in  agricultural  soils  in  Canada  is  generally  near  zero  (Smith  et  al,  1996;  Monreal  and 
Janzen,  1993;  Nyborg  at  al,  1995;  Liang  and  MacKenzie,  1992).  This  is  mostly  because 
changes  in  SOC  from  land  that  had  been  converted  to  agriculture  are  now  small  and  little 
new  land  is  being  converted.  Current  changes  in  SOC  are  mostly  a  result  of  changes  in 
management  practices  such  as  tillage,  crop  rotation,  and  amount  of  fertilizer  used. 
Management  practices  affect  SOC  storage  in  soils  by  influencing  the  amount  of  plant 
residues  returned  to  the  soil  and  the  rate  at  which  the  residues  decompose. 

The  work  carried  out  by  the  soil  degradation  indicator  group  has  been  aimed  at 
developing  methods  and  a  data  base  upon  which  trends  in  soil  degradation  processes 
may  be  tracked  through  time.   Indicators  for  soil  erosion,  salinity,  and  compaction  are  also 
being  developed.  The  information  from  these  studies  is  based  on  the  Soil  Landscape 
Database  of  Canada  (SLC)  and  the  1981  and  1991  census. 

Century,  a  site  specific  computer  simulation  model,  makes  use  of  simplified 
relationships  of  soil-plant-climate  interactions  to  characterize  the  dynamics  of  SOC  and 
nitrogen  for  grasslands,  crops,  forests,  and  savannas.   It  considers  several  agricultural 
management  practices  such  as  planting,  fertilizer  application,  tillage,  grazing  and  organic 
matter  addition.  Climate,  soil  and  erosion  data  are  used  in  the  model  as  input.  Century  is 
described  in  detail  in  Parton  et  al.  (1993,  1989,  and  1987).  It  has  been  extensively 
evaluated  under  contrasting  soil,  climate  and  agricultural  practices.  It  has  been  used  in 
Eastern  Canada  (Angers  etal.,  1993)  Western  Canada  (Monreal  era/.,  1995),  the  United 
States  (Parton  etal,  1987),  in  Northern  Europe  (Paustian  etal,  1992;  Parton  etal.,  1982), 
as  well  as  for  soils  in  the  tropics(Parton  et  al,  1989). 

The  objective  of  this  study  was  to  investigate  the  use  of  the  Century  model  for 
developing  a  SOC  indicator  for  Eastern  Canada.  The  Century  model  was  used  to 
estimate  the  rate  of  SOC  change  in  all  SLC  polygons  in  the  Mixedwood  Plains  Ecozone  of 
Ontario.    The  intensity  of  agriculture  in  the  region  and  the  large  amount  of  readily 
available  data  make  Southern  Ontario  an  excellent  location  for  developing  and  testing  a 
soil  organic  carbon  indicator. 


2.0  Methodology 

2.1  Modelling  Procedure 

The  Century  model  was  used  to  estimate  carbon  dynamics  for  262  SLC  polygons 
in  the  Mixedwood  Plains  Ecozone  of  Ontario  from  1960  to  2000.  Ten  polygons  in  the 
ecozone  were  not  included  because  of  missing  input  data.  The  runs  were  broken  into  two 
main  time  blocks  or  time  periods,  one  from  1960  to  1985  and  the  second  from  1986  to 
2000.  The  length  of  the  time  blocks  varied  slightly  to  accommodate  length  of  rotations. 

2.2  Model  Input 

The  Century  model  requires  the  input  of  soil  properties,  climate,  and  agricultural 
management  practices.  This  section  will  indicate  the  sources  of  input  data  and  describe 
how  the  data  were  used  for  the  model  runs. 

2.2.1  Soil  Properties 

Soil  properties  including  soil  texture,  pH,  bulk  density,  and  soil  organic  carbon 
were  obtained  from  a  recent  assessment  carried  out  in  Ontario  by  Bruce  Macdonald 
and  Bruce  Gleig,  of  the  Ontario  Land  Resources  Unit.  Their  survey  was  incomplete 
when  this  project  was  carried  out.  Missing  data  were  obtained  from  the  Soil 
Landscape  of  Canada  database.  Initial  soil  organic  carbon  (SOC)  was  determined  to 
the  30  cm  depth  from  the  A  and  B  horizons.  Fig.  1  shows  the  distribution  of  SOC 
within  polygons  of  the  Mixedwood  Plains  Ecozone  for  Ontario.  The  majority  of  the 
polygons  have  SOC  values  between  4.5  and  9.0  g  C  m"2. 

2.2.2  Climate 

Thirty-year  normal  monthly  minimum  temperature,  maximum  temperature  and 
precipitation  data  were  estimated  at  the  centroid  for  each  of  262  SLC  polygons  by 
"kriging"(a  technique  used  for  generating  maps)  the  weather  data  from  all  the  weather 
stations  in  and  surrounding  the  Mixedwood  Plains  Ecozone  in  Ontario. 

2.2.3  Crop  Rotations,  Fertilizer  and  Manure  Application 

Crop  rotations,  as  well  as  fertilizer  and  manure  application  rates  for  each 
ecodistrict,  were  determined  from  the  1981  and  1991  Census  (Bruce  Macdonald, 
personal  communication).  The  same  fertilizer  and  manure  application  rates  as  well  as 
crop  rotations  were  used  for  all  SLC  polygons  in  an  ecodistrict.  Crop,  fertilizer,  and 
manure  application  practices  for  1981  were  used  throughout  the  first  Century  time 
block  (1960  to  1985),  likewise,  practices  for  1991  were  used  throughout  the  second 
time  block  (1986  to  2000).  The  crop  rotations  used  in  this  study  are  shown  in  Table  1. 
Horticultural  crops  and  tobacco  were  not  included  in  the  simulations  because  these 
crops  are  not  included  in  the  Century  model. 
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Table  1.    Crop  rotations  for  the  Mixedwood  Plains  Ecozone 


Rotation 
Code 

Rotation 

R2 

Corn,  Beans,  Winter  wheat 

R4 

Corn  and  Beans 

R6 

Corn  (x2),  Spring  cereal  (x2),  Hay  (x3-5) 

R7 

Spring  cereal  (x2),  Hay  (x3-5) 

R8 

Continuous  corn 

R9 

Continuous  beans 

R10 

Continuous  hay 

R11 

Improved  pasture 

R12 

Unimproved  pasture 

Due  to  time  constraints,  carbon  dynamics  could  not  be  simulated  for  all 
rotations.  The  three  or  four  most  important  rotations  in  each  SLC  polygon  were 
simulated,  accounting  for  more  than  75%  of  the  area  in  27  of  the  30  ecodistricts.    The 
scheduling  of  events  including  planting,  harvest,  and  tillage  were  done  in  accordance 
with  Smith  et  al,  1996. 

2.2.4  Tillage 

Century  simulations  were  carried  out  for  conventional  and  no-tillage  systems. 
No-tillage  management  was  added  only  in  the  second  time  block  (1986  to  2000)  and 
simulations  were  carried  out  for  an  SLC  polygon  only  when  the  no-tillage  area  was 
greater  than  5%  (1991  Census). 

2.2.5  Erosion 


The  annual  amount  of  soil  erosion  in  1981  and  1991  was  based  on  an  indicator 
study  of  soil  erosion  risk  (Greg  Wall,  personal  communication). 


2.3  Calculation  of  the  Rate  of  SOC  Change 

The  rate  of  SOC  change  for  the  years  1981,  1986,  1990,  and  1991  was 
determined  by  taking  the  slope  of  a  10  year  regression  centred  on  each  particular  year. 
For  instance,  the  regression  for  the  year  1991  was  done  over  the  time  period  from  1986 
until  1996.  Such  a  procedure  was  required  to  account  for  rotations  which  were  several 
years  long.  For  rotations  of  durations  longer  than  5  years  the  regression  was  carried 
out  for  twice  the  length  of  the  rotation. 


iC=%F,(£FcCe) 


Calculation  of  the  rate  of  SOC  change  for  each  SLC  polygon  was  done  by 
weighting  the  results  by  tillage  type  and  crop  rotation  as  follows:  where  A  C=  rate  of 
SOC  change  for  a  polygon,  T  =  number  of  tillage  practices,  C  =  number  of  crop 
rotations,  F,  =  fraction  of  area  covered  by  tillage  practice,  Fc  =  fraction  of  area  covered 
by  crop  rotation,  and  Cc  =  rate  of  SOC  change  for  the  crop  rotation. 


3.0  Results  and  Discussion 

Almost  1400  runs  of  the  Century  model  were  carried  out  for  various  climatic,  soil, 
and  management  conditions  in  order  to  investigate  the  potential  of  using  this  model  as 
a  tool  for  developing  a  carbon  indicator .  As  indicated  in  Table  2,  the  Century  model 
predicted  that  there  is  a  small  amount  of  SOC  loss  in  the  Mixedwood  Plains  Ecozone  of 
Ontario  but  that  the  amount  of  loss  is  declining  over  time.  The  values  reported  in  this 
table  are  small.  Such  rates  of  change  cannot  be  measured  to  a  significant  degree  in 
field  trials  at  the  present  time.  Recent  studies  have  indicated  that  net  C02  emissions 
from  agricultural  soils  in  Canada  are  generally  near  zero  (Smith  et  al,  1996;  Monreal 
and  Janzen,  1993;  Nyborg  at  al,  1995).     This  is  mostly  because  changes  in  SOC  from 
land  conversion  is  now  small.  Current  changes  in  SOC  are  mostly  a  result  of  changes 
in  management  practices  such  as  tillage,  crop  rotation,  and  amount  of  fertilizer  used. 
In  Ontario  the  introduction  of  no-till  agriculture,  particularly  in  1991,  resulted  in  a  more 
rapid  decrease  in  SOC  loss. 


Table  2.    Yearly  rate  of  SOC  change  in  the  Mixedwood  Plains  Ecozone  of  Ontario 

Year  Rate  (kg  C  ha1)  SOC  (Mt) 
1981                                                 -39.6  -0.3 

1986  -34.9  -0.27 

1990  -32.1  -0.24 

1991  -24.3  -0.18 


The  rate  of  SOC  loss  of -24.3  kg  C  ha"1  in  1991  compares  reasonably  well  with  the 
value  of  -14  kg  C  ha"1  as  determined  by  Smith  et  al,  1 997.    The  rate  of  SOC  change  in 
this  study  was  calculated  from  262  polygons  whereas  it  was  determined  with  only  26 
polygons  in  the  study  by  Smith  et  al,  1997.  Soil  properties,  climate,  and  management 
practices  in  this  study  were  determined  from  more  complete  sources  than  in  Smith  et  al, 
1997. 

The  rate  of  SOC  change  within  each  ecodistrict  for  1981,  1986,  1990,  and  1991  is 
shown  in  Table  3.  There  is  quite  a  large  variation  between  ecodistricts.  A  wide  range  of 
SOC  rates  occurs  throughout  the  Mixedwood  Plains  Ecodistrict.  There  are  few  spatial 
trends.  There  is,  however,  a  greater  amount  of  loss  in  the  Southern  tip  of  Ontario  where 
agriculture  is  very  intensive.  Ecodistricts  557,  565,  567,  568,  570,  571,  and  572  all  show  a 
loss  of  SOC  from  1 980  to  1 991 .  Maps  of  SOC  change  at  both  the  polygon  and  ecodistrict 
scales  for  the  years  1 981 ,  1 986,  1 990,  and  1 991  are  shown  in  Appendix  A.   In  1 981 , 
Century  predicted  that  the  southern  tip  of  Ontario  and  the  area  southeast  of  London  had 
the  greatest  SOC  loss.  Considerable  loss  of  carbon  also  occurred  in  some  polygons  near 
Ottawa  but  other  polygons  in  the  area  gained  carbon  such  that  the  larger  scale  ecodistrict 
map  for  1981  shows  little  loss  of  SOC  near  Ottawa.  As  the  maps  progress  from  1981 
through  to  1991  less  SOC  is  lost  or  more  SOC  is  gained  throughout  the  Mixedwood  plains 
ecozone.  The  polygon  map  for  1991  still  indicates  that  greater  than  400  kg  ha"1  y"1  of  SOC 
was  being  lost  in  some  polygons,  particularly  southeast  of  London  where  the  soils  still 
have  a  high  SOC  content. 

A  frequency  distribution  of  the  rate  of  SOC  change  for  all  262  sampled  polygons  is 
shown  in  Fig. 2.   In  1 991  there  were  more  polygons  closer  to  equilibrium  than  in  previous 
years.  There  were  less  polygons  losing  or  gaining  large  amounts  of  SOC. 
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Table  3.    Yearly  rate  of  SOC  change  (kg  C  ha1)  for  each  ecodistrict  of  the  Mixedwood  Plains  Ecozone 
in  Ontario 

Ecodistrict  1981  1986  1990  1991  Polygons  Area  (kHa) 
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Figure  2    Frequency  distribution  of  SOC  change  for  262 
polygons  in  the  Mixedwood  Plains  Ecozone  in  Ontario 
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The  rate  of  SOC  change  for  each  crop  rotation  is  shown  in  Table  4.   Rotations  with 
corn  and  beans  lost  the  greatest  amount  of  SOC.  As  expected,  rotations  with  hay  and 
pasture  stored  the  most  SOC.  As  previously  mentioned,  more  SOC  was  lost  in 
ecodistricts  near  the  Southern  tip  of  Ontario  where  agriculture  is  more  intensive. 


Table  4.  Rate  of  SOC  Change  (kg  C  ha'1)  for  Each  Crop  Rotation 


Rotation 

Rotation 

1981 

1991 

Number  of 

code 

polygons  in 

which  rotation 

occurs 

R2 

Corn,  Beans,  Winter  wheat 

-1.02 

-14.8 

105 

R4 

Corn  and  Beans 

-53.3 

-66 

68 

R6 

Corn  (x2),  Spring  cereal  (x2),  Hay  (x3-5) 

-30.9 

-20.1 

234 

R7 

Spring  cereal  (x2),  Hay  (x3-5) 

-16.2 

-37.6 

86 

R8 

Continuous  corn 

-32.9 

-2.6 

56 

R9 

Continuous  beans 

-27.2 

-22.3 

10 

R10 

Continuous  hay 

5.8 

23.2 

29 

R11 

Improved  pasture 

32.7 

19 

168 

R12 

Unimproved  pasture 

30.8 

29 

199 

Soil  texture  strongly  influences  SOC  dynamics  in  agricultural  soils.  The  output  from 
Century  indicated  that  finer  textured  soils  usually  gained  more  SOC.   In  1991,  the  rate  of 
SOC  change  in  fine,  medium,  and  coarse  textured  soils  was  24.7,  0.6,  and  -78.8  kg  C  ha1, 
respectively.  Fine  textured  soils  were  classed  as  any  soil  finer  in  texture  than  a  silty  loam 
and  coarse  textured  soils  were  those  coarser  in  texture  than  a  loam. 
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Appendix 


Maps  of  soil  organic  carbon  change  in  the  Mixedwood  Plains  Ecozone, 
mapped  at  the  polygon  and  ecodistrict  levels 


Carbon  change  in  southern  Ontario  at  the  SLC  level  -  1981 
Carbon  change  in  southern  Ontario  at  the  SLC  level  -  1986 
Carbon  change  in  southern  Ontario  at  the  SLC  level  - 1990 
Carbon  change  in  southern  Ontario  at  the  SLC  level  -  1991 
Carbon  change  in  southern  Ontario  at  the  ecodistrict  level  -1981 
Carbon  change  in  southern  Ontario  at  the  ecodistrict  level  -1986 
Carbon  change  in  southern  Ontario  at  the  ecodistrict  level  -1990 
Carbon  change  in  southern  Ontario  at  the  ecodistrict  level  -1991 
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Rate  of  Change  (kg  /  ha  yr) 

■  <-4D0 
-100  to -CIO 

-loo  to  loo 

100  to  400 

|>400 

I      I  Data  incomplete 
to  rati  the  model 

1  Water 

urban  areas 


& 


Greenhouse  &  Processing  Crops 
R«*#4rcll  Cent**      |  Harrow,  Ontario  I 


Map  prepared  by 

Ontario  Land  Resource  Unit 

Guelph,  Ontario 


1*1 


■6C(t«:U.TU«Af» 
MMfOOD  CMWA 


Carbon  change  in  southern  Ontario  at  the  ecodistrict  level  -  1981 


^h^^r^yfl 


Lonaw 


a 


Rate  of  Change  (kg  /ha  yr) 

I        I -100  to -400 
I        1-100  to  100 
]  100  to  400 
^■>400 

nwitn 

I      I  urn  an  areas 


& 


Greenhouse  &  Processing  Crops 
RMMrch  C#ntr«      \  Harrow,  Ontario  I 


Map  prepared  by 

Ontario  Land  Resource  Unit 

Guelph,  Ontario 


1*1 


«Ri:ilLtVtt/«> 
AGBr-rOOD  DWVA 


Carbon  change  in  southern  Ontario  at  the  ecodistrict  level  -  1986 


Rate  of  Change  (kg  /  ha  yr) 


Zl -100  to -400 
I        1-lQQtOlCO 
Z\  100  to  400 

SB  >  400 

1      [water 

I      I  urti  an  areas 


^0 


& 


Greenhouse  &  Processing  Crops 

RMMrch  C#lttr#      \  Harrow,  Ontario ) 


M  ap  prepared  by 

Ontario  Land  Resource  Unit 

Guelph,  Ontario 


■♦■ 


AERrfXMD  OW1ADA 


Carbon  change  in  southern  Ontario  at  the  eco district  level  - 1990 


^f  pig- 


Rate  of  Change  (kg  /  ha  yr) 

|<-co 

100  to -400 

ioo  to  ioo 

]  100  to  400 
■  >400 

water 
urb  an  areas 


^ 


Greenhouse  &  Processing  Crops 
RftSMTCh  C#ltttft      (  Harrow,  Ontario  I 


Map  prepared  by 

Ontario  Land  Resource  Unit 

Guelph,  Ontario 


l+l 


Carbon  change  in  southern  Ontario  at  the  ecodistrict  level-  1991 
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